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5 ! Jet quenching of quarks and gluons in central A + A collisions suppresses the high-p T 

meson production and exposes novel baryon dynamics that we attribute to (gluonic) baryon 
junctions. This mechanism predicts baryon enhancement in a finite moderate-pr window 
that decreases with increasing impact parameter. We also extend recent calculations of the 
transverse momentum behavior of p/7r + and p/n^ to other baryon species and show that a 
similar pattern is expected for the K~/A and K + /A ratios. Within the framework of the 
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1. INTRODUCTION 

One of the unexpected results reported by PHENIX and STAR during the first year RHIC 
run at y/s NN = 130 GeV was that in contrast to the strong 7r° quenching for 2 GeV < pr < 



model constant p/p and A/A ratios are found in central Au + Au collisions at RHIC energies 
^ ■ up to p T ~ 4 — 5 GeV. 

O 
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^ ■ r un at y/ h NN 

(D \ 5 GeV, the corresponding charged hadrons were found to be suppressed by only a factor 

~ 2 — 2.5. Even more surprisingly, the identified particle spectra analysis at PHENIX 
suggests that Rb(pt) = vl^~ "> p/n + > 1 for p T > 2 GeV. Thus, baryon and antibaryon 
production may in fact dominate the moderate- to high-j>r hadron flavor yields [p]]. 

These and other data point to novel baryon transport dynamics playing role in nucleus- 
nucleus reactions. An important indicator of this is the high valence proton rapidity density 
at midrapidity (y — 0). More recent results corroborate the non-perturbative baryon produc- 
tion hypothesis through equally abundant A and A production [p]]. It has also been observed 
that the mean transverse momentum (pt)b for various baryon and antibaryon species is 
approximately constant and deviates from the common hydrodynamic flow systematics of 
soft hadron production in A + A collisions. This provides strong motivation to explore alter- 
native physical mechanisms that may give insight into the anomalous anti+baryon behavior 
at RHIC. 

Identified particle analyses from the recent y/s NN = 200 GeV RHIC run find similar 
puzzling features of moderate-pr baryon spectra that have to be understood in the context 
of even stronger quenching of neutral pions || . 
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Figure 1. The p/ir" ratio versus pr for Figure 2. The p/p ratio versus pt for cen- 
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Npart (solid) versus Nu n (dashed) geometry 
is shown. Boosted thermal source compu- 
tation and ratio of fits to PHENIX data are 
shown for comparison. 



130 GeV. Com- 
parison between a PQCD calculation and a 
quenched+baryon junction scenario. Cen- 
tral 10% data from STAR. 



2. BARYON PRODUCTION AND TRANSPORT MECHANISM 

A topological non-perturbative baryon production and transport mechanism was origi- 
nally proposed by Veneziano and Rossi in elementary nucleon-nucleon collisions ||. It has 
recently been successfully generalized and implemented for the case of nucleus-nucleus col- 
lisions P|- Its phenomenological applications are currently based on Regge theory, where a 
Regge trajectory J = a(0) + a'(0)M 2 is specified by its intercept a(0) and slope a'(0). It 
has been argued that a!j(0) ~ 0.5 and ctj(0) ~ l/3a^(0). Regge theory gives exponential 
rapidity correlations, which in the presence of two sources (at ±Y" max ) lead to net baryon 
rapidity density in central A + A collisions of the form: 



dy " ' sinh(l - aj(0))V max ' 

It is evident from Eq.(|l]) that the net baryon distribution integrates to 2 A and in going to 
peripheral reactions scales as N part . At RHIC energies of y/s = 130(200) GeV, corresponding 
to V ma x = 4.8(5.4), in central reactions dN B ~ B /dy = d(p — p)/dy + d(n — n)/dy + d(A — 
A)/dy + • • • ~ 18(13.5). The relative contribution of each baryon species can be evaluated 
from isospin symmetry and strangeness conservation (via comparison to midrapidity kaon 
production). 

Hadronic transport in small-to-moderate pt is effectively controlled by the slope of the 
Regge trajectory This would suggest that the baryon/meson mean inverse slopes in a 
phenomenological py-exponential (~ e~ PT ^ T ) soft particle production model are related as 
T B :T M ~ y/3: y/2. Soft pion production, however, is largely dominated by resonance 
decays, where the cumulative effect from the random walk in px due to string braking is 
destroyed. This leads to the relation (p T ) n : (p T ) K : (p T ) B ~ 1 : (1 y/2) : V3 (220 MeV : 
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Figure 3. Same centrality classes as in 
Fig. 1 for A/K + versus px- Au + Au re- 
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Figure 4. Same as Fig. 2 for A/A in central 
Au + Au at yfs NN = 130 GeV. Minimum- 
bias data from PHENIX. 



275 MeV : 380 MeV) |Q. One also notes that in the limit of pair production dominated by 
junction-antijunction loops (which we consider here) the transverse momentum distribution 
of antibaryons closely resembles that of baryons (with (px)§ = (Pt)b)- 

The high-pr part of the hadron spectra is computed from leading order PQCD augmented 
by the effects of initial multiple scattering (Cronin), nuclear shadowing, and parton energy 
loss. The most important modification to the perturbative scheme comes from the final state 
gluon bremsstrahlung computed here via the Gyulassy-Levai-Vitev (GLV) formalism 0]. 
Its centrality dependence determines to a large extent (especially for pions) the centrality 
dependence of the baryon/meson ratios. Similar technique has been used to estimate the 
A-induced initial parton broadening ||. 

3. RESULTS AND CONCLUSIONS 

We have studied the interplay between a novel baryon production and transport mech- 
anism (in rapidity y and moderate p T ) and the suppression of particle spectra as a 
result of the medium induced non-Abelian energy loss of jets H in Au + Au reactions at 
\fs NN = 130 GeV. We find an enhanced baryon/meson ratio Rb(pt) ~ 1 that decreases 
with centrality in a finite pt window as illustrated in Figs. 1 and 3. For px > 5 GeV the 
ratios reduce below unity and follow the perturbative calculation. We have extended recent 
studies of the pr-differential p/n + , p/^~ ratios {§ (Fig. 1) to K~ /A, K + /A (Fig. 3). Similar 
enhancement pattern if found at moderate 2 GeV < px < 5 GeV. Its centrality dependence 
is predicted to be weak. Figs. 2 and 4 show that in central reactions the predicted p/p and 
A/A ratios are p?- independent up to px — 4 — 5 GeV. At larger transverse momenta they 
decrease and approach perturbative estimates 0. We propose that the non-perturbative 
baryon dynamics may extend to px — 5 GeV. 

The enhanced baryon and antibaryon production at moderate transverse momenta ac- 
counts for the difference in the suppression factor Raa{pt) of inclusive charged hadrons and 
neutral pions 0] (see Fig. 5). We have demonstrated that when the non-perturbative baryon 
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Figure 5. Suppression of neutral pions and 
inclusive charged hadrons relative to the bi- 
nary collision scaled p + p result. Central 
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Figure 6. Same as Fig. 5 extended to large 
transverse momenta. Note Raa{pt > 5 — 
6 GeV) w const. 
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contribution becomes small Raa{pt) is approximately equal for all hadron species. This is 
illustrated in Fig. 6 for p T > 5 GeV. For the first RHIC run at \/s NN = 130 GeV up to 
p T ~ 10 GeV the suppression factor is found to have a value Raa{pt) — 0.3 and is almost 
independent of the transverse momentum ||. 

We have performed PQCD calculations at the highest RHIC energy of \/s NN = 200 GeV 
that also predict a constant suppression factor Raa{pt)- The detailed Raa(pt) behavior is 
shown to be a consequence of the combination of Cronin effect, shadowing, energy loss, and 
the computable shape of the underlying jet distributions. 
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